Observation of multifractality in Anderson localization of ultrasound 
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We report the experimental observation of strong multifractality in wave functions below the 
Anderson localization transition in open three-dimensional elastic networks. Our results confirm 
the recently predicted symmetry of the multifractal exponents. We have discovered that the result 
of multifractal analysis of real data depends on the excitation scheme used in the experiment. 
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Critical phenomena are of prominent importance in 
condensed-matter physics. Criticality at the Anderson 
localization transition has been the subject of intensive 
theoretical research. Some important theoretical predic- 
tions have been made, among which is the remarkable 
aspect of multifractality of wave functions at this tran- 
sition. Numerical simulations support these predictions 
but also raise more questions Recent experimental 
progress has paved the way for the direct investigation of 
the Anderson localization transition at the mobility edge 
in real samples @, S 0] ■ 

In this Letter, we report the experimental observation 
of strong multifractality (MF) just below the Anderson 
transition. This observation is based on the excitation 
of elastic waves in an open 3D disordered medium. The 
recently predicted symmetry relation of the MF expo- 
nents [fj is tested and confirmed. All results are com- 
pared with the corresponding analysis of diffusive (metal- 
lic) wave functions in the same network at a different 
frequency or with a light speckle pattern generated by a 
strongly scattering medium, showing a very clear differ- 
ence between localizing and diffusive regimes. Our results 
not only highlight the presence of MF in wave functions 
close to the mobility edge, but also reveal new aspects of 
the MF character in real experimental systems. 

Before presenting the experimental results, we briefly 
review some general aspects of MF and their implications 
in the context of the Anderson transition. Multifractal- 
ity quantifies the strong fluctuations of the wave func- 
tion. It shows the non-trivial length-scale dependence of 
the moments of the intensity distribution. The depen- 
dence can be investigated by varying the system size L, 
or alternatively, if the system size is fixed, by dividing 
the system into small boxes of linear size b and varying 
b. This property is quantified by using the generalized 
Inverse Participation Ratios (gIPR) 
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where I(r) is the normalized intensity (equal to 



IV' 2 ( r )l/ / \ip 2 (^)\d d r where ip(r) is the wave function) 
and Jb 4 is the integrated probability inside a box Bi of 
linear size b, with A -C b <C L where A is the wavelength. 
The summation is performed on the whole sample, which 
consists of n — (L/b) d boxes, and d is the space dimen- 
sion. By definition P± = 1 and Pq = n. 

At criticality, the ensemble averaged gIPR, (Pq), scales 
anomalously with the dimensionless scaling length L/b as 
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where d(q — 1) and A g are called the normal (Euclidean) 
and the anomalous dimensions, respectively. For a nor- 
mal (extended) wave function, A g = for every q. A 
(single-) fractal wave function with fractal dimension D 
is described by r(q) = D(q — 1). For critical states r(q) 
is a continuous function of q that fully describes the MF. 

MF describes the scaling of the moments of a probabil- 
ity density function (PDF). The gIPR, defined in Eq. (TT|), 
are proportional to the moments of the distribution func- 
tion of the eigenfunction intensities, so that Eq. ((2|) im- 
plies the following scaling relation for the PDF: 
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The second term in the exponent, f(a), is called the sin- 
gularity spectrum, and is related to the set of anomalous 
exponents r(q) by a Legendre transform 

r(q) =qa-f(a), q = f'(a), a = T '(q). (4) 

The singularity spectrum f(a) is the fractal dimension 
of the set of those points r where the wave-function in- 
tensity, I(r), scales as L~ a . In mathematical terms, it 
shows the coexistence of several populations of singular- 
ities in the measure, which is the wave-function inten- 
sity for this specific case. In the field-theoretical treat- 
ment of random-Schrodinger Hamiltonians, MF implies 
the presence of infinitely many relevant operators [y, 0] • 
The functional dependence of f(a) is an important and 
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unique property of each universality class. In the ex- 
tended regime, V(IiiIb) is strongly peaked near a = d, 
since the short-range "Gaussian" fluctuations [|| are 
washed out in the box integration. 

First order perturbation theory for an Anderson transi- 
tion in 2+e dimensions (two is the critical dimension), 
results in the "parabolic approximation" for the MF wave 
functions [l(J [llj] . This result, A g = 75(1 — q), corre- 
sponds to /(a) = d — (a — d — 7) 2 /Aj, where 7 is a con- 
stant in the order of e. A similar approximation app lies 
to metallic (diffusive) states in three dimensions [12 ll3T ] 
due to weak localization, although with 7 < 1. This is 
sometimes called weak MF. 

Recently, an exact symmetry relation 
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was theoretically predicted for the set of anomalous ex- 
ponents [f|. The numerical and analytical investigations 
of the 3D Anderson model and certain random matrices 
suggest that MF may exhibit itself also for off-criticial 
states on both sides of the transition 14] . The MF con- 
cept was extended to the boundaries where it behaves 
differently with respect to the bulk (l5j . 

Most of the available information about MF is based 
on numerical investigations (See, e. g., [l|, [l|| 12, 0| 
and references therein). The only experimental attempt 
to observe strong MF in wave functions so far is due to 
Morgenstern et. al using scanning tunneling microscopy 
of 2D electron systems [2| . Their observation of MF was 
hindered by the presence of several eigenfunctions in the 
measurement and by the limited size of their system. 

We now use ultrasonic measurements to demonstrate 
three different, but closely related, manifestations of MF: 
(1) the probability density function (PDF), (2) the scal- 
ing of generalized Inverse Participation Ratios (gIPR), 
and (3) direct extraction of the singularity spectrum. 
Our experiments were performed on disordered single- 
component elastic networks, made by brazing randomly- 
packed 4.11-mm-diameter aluminium beads together [3j. 
The data presented here were obtained from a represen- 
tative disc-shaped sample with a 120 mm diameter and 
14.5 mm thickness. Two different configurations were 
used for excitation. In the first excitation scheme a point- 
like ultrasound source emits short pulses next to the sam- 
ple surface. In the second case the source was put far 
from the sample so that a quasi-planar wave was inci- 
dent on the whole interface. In close proximity to the 
opposite interface, vibrational excitations of the network 
were probed with sub-wavelength-diameter detectors in 
the frequency range of 0.2 to 3 MHz, where the wave- 
lengths are comparable to the bead size and the scatter- 
ing mean free paths are much less than the sample thick- 
ness [3j. The intensity at a particular frequency was de- 
termined from the square of the magnitude of the Fourier 
transform of the entire time-dependent transmitted field 
in each near-field speckle. The intensity was normalized 
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FIG. 1: Comparison between coarse-grained PDF for local- 
ized and diffusive speckle intensities. The PDFs are exper- 
imentally obtained from the histogram of the logarithm of 
averaged intensities in the localized (thick bars) and diffusive 
(thin bars) regimes. The black line shows a fit to a single 
parameter log-normal distribution given by the parabolic ap- 
proximation of Eq. [3] Inset: The peak position (symbols) 
and the full width at half maximum (bars) of the intensity his- 
togram is plotted for localized (circles) and diffusive (squares) 
speckles as a function of coarse-graining box size. 



by the total intensity in the measured speckle pattern. 
The normalized speckle intensity, I (J) was recorded at 
each point j on a square grid of linear size L g = 55 points 
with a typical nearest-neighbor spacing of 0.66 mm. 

In the lower frequency band around 250 kHz, the ul- 
trasound propagation is diffusive. A localizing regime is 
observed in a 50% bandwidth around 2.4 MHz, where 
the measured localization length in the sample is smaller 
than the size of the analyzed speckle patterns (0.7L S ) 
and almost equal to the sample thickness. A full de- 
scription of the experiment and a thorough comparison 
of previous measurements with the self-consistent theory 
of localization has been presented in Q . 

We obtain the PDF from the histogram of the loga- 
rithm of box-integrated intensities Jb 4 . We sample over 
100 speckles in a 5% bandwidth around 250 kHz and 
2.4 MHz for diffusive and localized regimes, respectively. 
Two representative histograms are shown in Fig. [1] with 
typical box sizes of b = 9 and 6 = 2 points for low and 
high frequency measurements, corresponding to box sizes 
of approximately two wavelengths in both cases. The 
PDF for localized waves is clearly much wider than the 
one for diffusive waves and the peak is shifted from the 
average intensity. We have also plotted the peak-position 
and the width of the histogram as a function of box size 
in the inset of Fig. [TJ 

In principle, it is possible to extract the MF spectrum 
from the PDF [l8[. However, a box-counting analysis 
can give more accurate results based on the scaling of 
the gIPR. Similar to many numerical studies, we approx- 
imate the expectation values by box-sampling over a sin- 
gle or multiple wave functions measured for a single re- 
alization of disorder. This approximation is known as 
typical averaging. Typical averaging is unable to reveal 
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FIG. 2: (a) The measured anomalous exponents A 9 are shown 
for localized ultrasound (full squares) and diffusive light (open 
circles) speckles. The dashed line shows the same data-points, 
mirrored relative to q — ~ in order to check the symmetry in 
the spectrum. The anomalous exponents are estimated from 
the box-counting method. The slope of the gIPR plotted ver- 
sus the box size in bilogarithmic scales yields A,. One ex- 
ample is shown in the inset for q G {—2,-1,0,1,2,3} and 
/ = 2.40 MHz. (b) The average singularity spectrum is cal- 
culated for the ultrasound speckles (full squares) at frequen- 
cies between 2.0 to 2.6 MHz. For comparison a singularity 
spectrum for diffusive optical speckle (open circles), with the 
Euclidian dimension, is extracted by applying the same box- 
counting procedure. 

information that is related to statistically rare events [l| . 
In this approach, the system size is fixed and supposed 
to be large enough relative to the box size. The approxi- 
mate scaling relation is derived by plotting the estimated 
gIPR, given by Eq. (J]), versus the box size b [24(. Note 
that although we have examined three-dimensional sam- 
ples, the Euclidian dimension of our sampling space is 
two since the available data is taken just from the sur- 
face of the sample. The effective system size is L g over 
which the intensities are normalized. By plotting P q ver- 
sus the box size in bilogarithmic scales [e.g., see the inset 
to Fig.^a)], power-law behavior is found for q <G [—3,4], 
with the slope yielding the scaling exponent r(q). The 
average anomalous exponent is obtained by averaging the 
exponents measured for several frequencies between 2.0 
and 2.6 MHz and subtracting off the normal part of the 
exponent 2(g — 1). The standard deviation is taken as 
the error-bar. 

The anomalous exponents are plotted as a function of 
q in Fig. Ufa). For comparison with the localized data, 
the same numerical procedure was applied to a diffusive 
speckle pattern, where the behavior is entirely different 
(A q = 0). In making this comparison, an optical diffu- 
sive speckle pattern was used to capitalize on the best 
available statistics. 

The behavior of the anomalous exponents shown in 
Fig. [5] provides unambiguous evidence for surface MF of 
the localized ultrasound wave functions. This is the most 
important result in this Letter. Note that MF is clearly 
seen in these data, even though the localized wave func- 
tions in our finite sample are near to, but not exactly 
at, criticality. In addition, our observation of MF clearly 
supports the predicted symmetry relation ([5]). Our ex- 
perimental demonstration of this fundamental symmetry, 




FIG. 3: (a) The reduced anomalous dimension 8 q = ^7rj^) 
is plotted versus q. Bars show the estimated error. Devia- 
tion from a horizontal line corresponds to the deviation from 
parabolic approximation, (b) The reduced anomalous dimen- 
sion 5i is plotted versus frequency in the localized regime 
for two excitation schemes: point-source (squares) and plane 
wave (circles). The error bars represent the standard devia- 
tion of the measured exponents that are averaged over each 
0.1-MHz- wide frequency band. 

seen in a very different system to the ones envisaged in Q , 
attests to the universality of critical properties near the 
Anderson transition. 

Finally, we have extracted the surface MF spectrum 
directly from the measurements fl9| . In this method the 
numerical error caused by the Legendre transform (j4]) is 
avoided. To get enough statistics, 100 wave functions in 
a bandwidth of 5% are used to estimate the MF spec- 
trum for several seven frequency bands between 2.0 and 
2.6 MHz. No systematic deviation is observed between 
the seven spectra obtained in this frequency range. These 
spectra are then averaged for each value of q G [—6,6] and 
the standard deviation is considered as the error bar. The 
results are summarized in Fig. [2b). The peak of the MF 
spectrum is shifted from two (the Euclidian dimension of 
the measurement basis) by a value of 0.21 ± 0.02. For 
comparison, the same procedure is applied to the optical 
speckle using the same g-range. No shift is observed for 
the optical speckle. 

The MF that is clearly seen in our data allows us 
to test the deviation from the parabolic approximation. 
This is characterized by the reduced anomalous expo- 
nents S q = q(il q ) • I 11 our results, shown in Fig. GHa), 
we see a deviation of less than 20% for q £ [—3,4]. The 
non-parabolicity of the spectrum is very difficult to mea- 
sure but it may have important theoretical consequences. 
More precise investigation of larger samples is needed to 
reliably confirm or exclude the possibility of a small but 
significant deviation. 

We have also investigated the dependence of the re- 
duced anomalous exponent at the symmetry axis, 5i — 
4Ai , on the frequency and type of excitation. The results 
are presented in Fig. [3]Jb). We observe a robust pres- 
ence of MF for all frequencies between 1.7 to 2.9 MHz. 
The measured anomalous exponent is larger for the point 
source illumination. This difference may be related to 
the number of modes excited in each scheme. It has 
been previously discussed [2(| that the overlap of two or 
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more eigenmodes shifts the peak of the singularity spec- 
trum towards the Euclidian dimension. Since the surface 
area of the sample is larger than the localization length, 
neighboring localized modes may coexist at the same fre- 
quency. These modes can all be excited by a quasi-plane 
wave while a point source couples more efficiently to the 
closest mode. 

Numerical analyses of bulk and surface MF for the 
eigenstates of the Anderson tight-binding Hamiltonian 
on a 3D cubic lattice at the metal-insulator transition 
have predicted corresponding shifts of 1.0 and 1.6 from 
the Euclidean dimension for the peak of the singularity 
spectrum 18|, l2l| . Another numerical study for an equiv- 
alent vibrational model on the fee lattice shows a similar 
outcome for bulk MF 22j, indicating to the universality 
of this phenomenon. It is not simple to explain the dif- 
ference between the available numerical results and our 
experimental outcome. Several issues may play a role. 
Mode overlap and the finite lifetime of modes due to open 
boundaries are two of these issues. Most numerical stud- 
ies are done based on uncorrelated disorder, which is ex- 
perimentally hardly ever achieved. The uniform bead size 
in our samples, which is comparable with the vibrational 
wavelength, is a source of correlation. The presence of 
correlation in the disordered potential may influence the 
critical behavior and induce nonuniversality [23I ]. 

Despite the wealth of theoretical and numerical stud- 
ies on the Anderson transition in 2D and 3D for the 
Schrodinger Hamiltonian in closed systems, critical prop- 
erties of this transition for classical waves in an open 
system have never been studied. Our system is specially 
challenging due to its 3D nature, open boundaries and 
presence of three polarizations for the elastic waves. Spe- 
cific properties of classical waves such as absorption are 
yet to be investigated in the context of criticality. Our 
results show that these important questions can now be 
investigated experimentally, providing vital guidance for 
new theoretical work. Our experiments reveal that the 
concept of MF not only concerns critical states but is 
valid as well around the mobility edge. This observation 
agrees with recent theoretical investigations [3] • Mutual 
avoidance of wave functions at large energy separations 
and their enhanced overlap at small energy separations 
are other important predictions of that theory, which can 
also be verified by our experiment. 

In conclusion, we have presented the first experimen- 
tal observation of multifractal wave functions below the 
localization transition. Our data validate experimentally 
the predicted symmetry relation of the anomalous expo- 
nents. Free from interactions and with the possibility of 
diverse illumination and detection schemes, sound and 
light experiments can provide a tremendous amount of 
useful information for this field of research. We believe 
that our observation of multifractality in classical waves 
will stimulate new theoretical and numerical investiga- 
tions. On the experimental side, this work highlights 



again the strength of statistical methods for studying An- 
derson localization. 
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